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ABSTRACT 


A lifting-surface theory and numerical procedure for 
designing supercavitating propellers are presented. 

Both a subcavitating and a supercavitating propeller 
are represented by vortex and source distributions. Unlike 
the subcavitating propeller, however, source strengths for 
a supercavitating propeller are related to cavity thick- 
ness, which is not known without examination, and have to 
be obtained by solving related integral equations. 
Numerical solution of the integral equations is obtained as 
a correctional function for a stripwise supercavitating 
cascade theory which, with lifting-line theory, is used for 
preliminary design of supercavitating propellers. The in- 
duced axial, radial, and tangential velocities are obtained 
on a blade reference surface that allows arbitrary skew, 
rake, and radial pitch variations. The blade shape is 
obtained as a correction to the shape obtained from strip- 
wise supercavitating cascade theory. Thrust and torque 
coefficients are obtained from pressures on the blade 
surface. 

The method is applied in designing several super- 
cavitating propellers that have design conditions the same 
as those of existing supercavitating propellers. 

Numerical computations were also used to design two 
additional supercavitating propellers which were built and 
tested. The design predictions are compared to the experi- 
mental data, both for blade cavity height and performance 
characteristics, and good correlation is obtained. 


ADMINISTRATIVE INFORMATION 

The Naval Material Command (NAVMAT O08T) funded this investigation in support 
of an ongoing Ship Performance and Hydromechanics Exploratory Development Program 
(Program Element 62543N, Task Area ZF43421001) assigned to the David W. Taylor Naval 
Ship Research and Development Center (the Center). This work has been completed 
under the High-Speed Propulsor Task portion of this program, (Center Work Unit 
1500-103), and partly under the Naval Sea Systems Command's General Hydrodynamic 
Research Program (Center Work Unit 1542-817). 


INTRODUCTION 


Ever since supercavitating propellers were first investigated systematical- 
1 


‘3 


ly, steady progress has been made in understanding associated problems. 


Designers of supercavitating propellers have been urged to consider the interference 


*A complete listing of references is given on page 83. 


effect of neighboring cavities and the effect of flow recacdations sou To this end, 
three-dimensional cavity-flow theory similar to the lifting-surface theory for sub- 


7-10 


208 has been formulated. However, to date none of the - 


cavitating propellers 
theories has been applied to the actual design of a supercavitating propeller owing 
to the numerical complexity of the problem. 

As with subcavitating propellers,” a supercavitating propeller may be designed 
in two steps: preliminary or lifting-line design and final or lifting-surface 
design. The former gives an approximate solution and supplies basic data to the 
latter. The preliminary design supplies information for the final design concerning 
load distribution on the blade and preliminary pitch distribution, which forms the 
basic singularity surface for the final design. The section shape of the blade 
cavity and the cavity length of a supercavitating cascade model are also computed in 
the preliminary design. The information is used in the final design as the first 
approximation. Although this report deals mainly with the final design, we also 
consider aspects of the preliminary deeienn to help the reader understand the final 
design more fully. 

The propeller diameter, the blade number and contour, the hub diameter for the 
given propeller-thrust power coefficient with the given ship speed, propeller revo- 
lutions per minute (rpm), and the wake fraction are considered to have been deter- 
mined before applying the present preliminary and final programs. To ensure that 
the designed propeller is fully cavitating and that the blades are thick enough, the 
leading-edge cavity thickness is specified. If the specified cavity thickness near 
the leading edge does not accommodate a stable cavity of at least 1.5 chord lengths 
because the cavitation number is too large, then the leading-edge cavity thickness 
is increased internally to make a long enough cavity for the preliminary program. 

The purpose of the present program is to design a fully cavitating propeller 
that is efficient, has no face cavitation, is structurally strong, and meets the 
design requirements. 

As for a subcavitating propeller, the supercavitating propeller is represented 
by vortex and source diseribucioast. The source strengths are related to the cavity 
thickness, are not known a priori and are to be obtained by solving related singular 
integral equations. When the source strengths are known, the supercavitating pro- 
peller problem is similar to the problem for a subcavitating propeller with thick 


and wide blades. Thus, many of the computation techniques apply to both 


5; 


subcavitating and supercavitating propellers. Indeed, this present study makes 
use of many parts of the lifting-surface design programs reported by Kerwin for 

a subcavitating propeller. The coordinate systems of the blade are the same for 
both propeller types. 

Since the singular integral equation is a Fredholm equation of the first kind, 
the method for its solution must be chosen with extreme care. Additionally, the 
cavity shape is not known without examination. Thus, iteration and/or some special 
cavity model has to be applied knowing that all the inviscid cavity models are not 
exact; rather, they are approximate representations. The source distribution is 
obtained from a stripwise, two-dimensional, supercavitating-cascade representation 
developed in the preliminary design. For the procedure used at present, the source 
distribution is multiplied by a double polynomial having unknown coefficients in 
terms of the radial and chordwise coordinates. Since the two-dimensional cavity 
model used in preliminary design is a linear, double spiral vortex model, the cavity 
streamline is closed at infinity instead of at the cavity end. Thus the unknown 
source strengths are distributed in a prefixed plane that contains the blade surface 
and extends to the wake about 1-1/2 chord lengths. At the end of the plane is 
attached another source line having an unknown polynomial strength that is to be 
solved along with the double polynomial coefficients by the least squares method. 
The solution is a function of propeller geometry, given thrust or power, blade load 
distribution, advance coefficient, and cavitation number. 

The induced axial, radial, and tangential velocities--thus pitch and camber 
distribution--are obtained on a blade-reference surface that allows arbitrary skew, 
rake, and radial pitch variation. The blade cavity shape is obtained as a correction 
to the blade cavity derived from supercavitating cascade theory. The thrust and 
torque coefficients are obtained from the pressure on the blade surface, which is 
converted from the lift distribution through the Kutta-Joukowsky theorem. 

If the hub boundary condition is not considered in solving for the cavity 
source distribution, the radial velocity caused by the cavity source is not stable. 
When the hub boundary condition is considered, the numerical results are shown to 
have reasonable convergence in many numerical experiments. Two new supercavitating 
propellers were designed using numerical computations derived during the present 


program, and models were manufactured. Model test results correlate reasonably well 


with the computed cavity shapes and powering performance. 


GEOMETRY OF THE BLADE 

The geometry of the blade of a supercavitating propeller is almost the same as 
that of a subcavitating propeller. The blade is represented in the flow by singu- 
larity distributions of vortices and sources on the blade-reference surface in the 
linear version. Thus, the blade-reference surface has to be known approximately, 
although the blade shape has to be obtained as a solution. As in the theory for sub- 
cavitating propellers, the reference surface will be close to the helical surface of 
the hydrodynamic pitch angle Be obtained in the preliminary design. Although the 
reference surface could be determined more accurately by an iterative procedure, in 
the present program Be is used once without iteration, which is considered to be a 
reasonable approximation. The pitch angle of the reference surface of the wake can 
be set differently from that of the blade. 

A cylindrical coordinate system (x,r,8), seen in Figure 1, is defined; the x- 
axis is coincident with the axis of rotation of the propeller and is positive when 


facing downstream. Thus, the reference surface can be represented by 


xm="ANOi+ Gr), ty See at ails oF <0<o (1) 
\ =F tan’ Br (2) 
al 


where T(r) depends on rake and skew, and 4 is in general close to a constant. If a 


Cartesian coordinate is defined as in Figure 1, the parametric representation of the 


reference surface x(r,9) = xi + yj + zk is 
x(G- 6), = 9G tear) 
y(c,0)) = Eacos 0 (3) 
z(r,8) = r sin 6 


Then the area element is 


ds = /EG-F* drd® = H drdé (4) 


where 
E = Xp ° Xp = ae ees Be 
ete eee Peer een aaa (5) 
E x. r es 
F = Se i A(8A +E) 
Thus 
1/2 
as = [+e +e? (Ate) dae (6) 


The unit vector normal to the reference surface is 


x. x x. : Xp x x ; (-ritr (A 8+0,)I+\k os 


10 lee H 


= = 

xp, | EG-F* 
If X} and © are independent of r, the normal vector in Equation (7) is on the cyl- 
inder, r = constant. That is, the radial component of the normal vector is zero. 


Two convenient unit vectors tangential to the reference surface are 


ete) 
Sm (8) 


and 


(9) 


The curvilinear distance is measured along the helix on the reference surface 
from the 6 = O plane to the leading and the trailing edges, denoted by So and Sie 
respectively. The chord length of the propeller blade at any radius is 


c(r) = S(r) = S, (r) (10) 


VELOCITIES DUE TO SINGULARITIES 
Both a source distribution with strength m(r,8) per unit area and a vortex 
distribution with strength G(r,6) per unit area are distributed on each blade- 
reference surface. Then, the velocity u" caused by the source distribution is 


derived from potential theory 
—m yal (=) 
u (x,r;0) aX At {aco ie B H dod (it) 


The velocity ae caused by the vortex distribution is derived from the Biot- 


Savart law 


a (a,2,0) = fe | OC.) SES (12) 


where 


Die ten 2) 1/2 
(x-§) +r +o -2rp cos (b45,-0)| (13) 
o 2e0e1) 
é. = 7 (14) 


B = {x-E(p)}i + {r-p cos (46, -8)}5 - {p sin (o-8,-8) }k (15) 
where Z is the number of blades 
d2 = {tan Si+ cos (+5, -8) 5+ sin ($+6,-9)k}dp (16) 


which is a line element along the intersections of the 9 = ¢ plane and the reference 


surface. 6 is the angle between d& and the x = constant plane or tan 6 = 0x/dr = Or. 


+ (4). The axial, tangential, and radial components, us> Ug> and u are written as 
follows 
iL oo i 
m Ak 1 
u (x, r0)) = pes m(p,o) aay (x=€) H dodp 
a 4 D) Bo 
Ty oF k=1 
1 o@ Z 
Te Gene tae { m(p,¢) y 2 sin ($+6,-6) H dodo 
Oi a? At i 3 k 
8 k=1 ® 
EGS a 
1 © Z 
m ut iE 
u_(x,r,6) = =— i m(p,) s — {r-p cos ($+6,-6)} H dodp (Glyp) 
r At 3} k 
9 ki 8 
Pilea 


al on Z 
bene 0), = é | i G(p,¢) s + {x sin (p+5,-0)]} H dbdp 
a TT B 

Ty a k=1 


(cont. ) 


IL 00 i, 
: are, 
ss All I G(0,9) s eS) {p"-pr cos ($+6,-6)} H dodo (18) 
ry OL k=1 


ts) 


AV Vb 
Gc ale ilk 
esa) So | i CO oe es en ee) 
: ii 78 pa 


H L 


x cos (+6, -6) }-p tan 6 cos (p45, -6)] H dodo 


1 oo Z, 
A i | i C.(0,9) > “3 [p tan B, cos (p+5,—-8) 
a oF = 
2» (x-&) sin ($46, -8)-r tan B] H dodo (18) 


il On ZL 
ue (x41) => { | G(p,%) y +, {(x-&) sin (o+6,-8) 
ly k=1 : 


a 
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where Ty = hub radius 
Ce = - dG(p,6)/do = trailing vortex strength in oF <o¢< © 
subscripts a, 8, and r indicate axial, tangential, and radial components, 
respectively 


The induced velocity component due to the singularity distributions, normal and 
tangential to the blade-reference surface, can then be written in a linear approxi- 


mation from Equations (5) and (7) 


ee if i ie dA C 
us os HG,6) ru, A(r)uy T (oe 4 + as) u. (19) 
ages eS en) (20 
= NyHGye) vat 6 ) 
where 
Hy =r sec BS (20-1) 
CAVITY BOUNDARY CONDITION 
The cavitation number defined with respect to ship speed Nie is written 
Jee let 
00 e 
Oo = (21) 
le bee 
2 p s 


where P is the pressure infinitely far upstream and Po is the pressure on the 
cavity. From the Kutta-Joukowsky theorem the vortex distribution can be written in 


a nondimensional form (Appendix A) as 


72 (p+0) Zt (22) 
where 
: P - Ee 
24262) (23) 


P is the pressure on the foil, and Vo is the local velocity. In designing a foil, 


we can use the local velocity obtained during preliminary design 


2 2 1/2 
My {(V+u,) +(rutu,) } (24) 


where V = V_ (l1-w) 
a s 
w = wake fraction 
® = angular speed 


When the supercavitating propeller, advancing with velocity ves is represented 
by source and vortex distributions, the perturbation velocity component parallel to 
the blade-reference surface can be written by a linear approximation of the Bernoulli 


equation of flow, referring to a moving frame of reference (Appendix A), as 


m G 
Un Ur pes 
Ver yore 2 (2) 
s s 
on the pressure side of the blade, and 
m G 
u u Vv 
ag, oil 2,(6h 8 
Venys a008G eo 
Ss s 
on the cavity, since 
G Gree 
Uny = Up =G 
(27) 
G G2 
Uy af Un 2 Ur 


where Up is the tangential velocity due to the vortex distribution at all points 


except the point concerned; subscripts "plus" and "minus" indicate the values on the 
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upper and lower sides of the mean camber of the blade, respectively; from Equations 
(25) through (27) the tangential velocity due to the source distribution alone can 


be written as 


We Vv V u 
i G a Oiegs f£ 
ye a Daye + Tew Ue on the blade plane (28) 
Ss s s 
and 
un Oo ie Ue 
Va oD Wels on the cavity plane (29) 
s s behind the blade 


When G and o are known, these two equations will form a system of Fredholm integral 
equations of the first kind for the source distribution. 

When the propeller blade shape is given instead of the vortex distribution, the 
problem becomes one of prediction rather than design, and the boundary condition on 


the pressure side of the blade is 


ee (30) 


where F is the blade ordinate with respect to the reference surface and x = pg sec 
Bo. 

a 

In this case the blade reference surface must be in close agreement with the mean 
camber and the wake surfaces. Now Equations (28) through (30) become the integral 
equations for both sources and vortices. Thus, the prediction problem becomes more 


complicated. In this report only the design problem is considered. However, the 


prediction problem can be solved similarly. 


SOLUTION FOR SOURCE 
For the lifting-surface design of supercavitating propellers, the load distri- 


bution p +o on the blade is supplied from the preliminary design computations. As 


a 


for subcavitating propellers, the pitch angle B of the blade-reference surface is 
taken as the hydrodynamic advance angle given by the preliminary design calculations. 
Then, were the source distributions known, all the induced velocities could be 
computed, and the blade-cavity shape and the final pitch could be obtained in a 
manner similar to that used for a subcavitating propeller. 

Two problems, strength and the location, are related to the source distribu- 
tions. The location may be considered to be on the blade reference surface. But, 
how far downstream should the sources be distributed? A two-dimensional supercavi- 
tating cascade section could be considered to contribute to the source distribution. 
When the Riabouchinsky cavity model is used, the cavity source will be confined 
inside a finite cavity domain. Thus the cavity domain supplied by two-dimensional 
theory could be used in an attempt to solve the source strength distribution aeiee 
If the linear double spiral vortex model is used, the source distribution does not 
terminate at the cavity end, but the wake source continues from the end of the 
cavity. Although there exists a logarithmic singularity in the normal velocity at 
the cavity end, the streamline represented by the integration of the normal velocity 
continues smoothly at the cavity end. Any cavity mode is known to give reasonable 
predictions of lift and drag. Therefore, it may be better to free ourselves from the 
concept of a finite cavity platform and to consider the cavity wake-source distribu- 
tion as beginning behind the leading edge. The major question is where to truncate 
and how to reduce the truncation error. 

Even if the problem of the source domain were solved, the method of solution 
would be far from definite because there is no established method to solve the 
Fredholm singular integral equation of the first kind; there are suggestions, how- 
ever, that a series of eigenfunctions could be sede Thus solutions to the three- 
dimensional problem of a supercavitating foil are complicated and time consuming. 

To arrive at a reasonable method to solve the present problem, we considered a 
combination of three-dimensional corrections to a two-dimensional solution and the 
series of eigenfunctions. That is, the two-dimensional solution for the supercavi- 
tating cascade is multiplied by a double polynomial function with unknown coeffi- 


cients 


J J adh 
£(p,$) = ) ai eee at > S aa Oa sae in the blade plane (31) 
ail i=l j=2 


with x = Ho, which is zero at the trailing edge, and 


J) al) J 
S i-l s ’ S ’ cee teat Dn 
(01,0) = 454 oO” an bay onl x" in the cavity plane (32) 
i=l i=l j=2 


The two-dimensional solution is supplied from preliminary design of each 
section. The cavity plane has to be finite numerically. A separate source distri- 
bution along the truncation line of the cavity is considered with unknown coeffi- 


cients to compensate for truncation of the cavity: 


i} 
f,(0) = wag Da, o 7 (33) 
i=1 


In addition, the cavity thickness near the leading edge is specified to be equal to 
that computed from cascade theory during preliminary design. This thickness is 
chosen because many sets of cavity thicknesses give the same load distribution and 
the same cavitation number for the linearized, two-dimensional cavity problem. The 
leading edge conditions can be interpreted as follows: the source distribution at 
the leading edge is exactly the same as that of a two-dimensional supercavitating 


cascade, that is 


Ee Pe the leading edge ee ee 


With the previously described representations of the source distribution, the 
integration of each term in the integral equation is performed numerically. If the 


propeller blades have rake or skew, & is a function of p. This is treated by a 


linear approximation in each Ap interval 


E = BG © f (34-1) 


where 
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If } = p tan B. is assumed to be a constant in interval Ap, all the integration with 
respect to p can be performed analytically in the interval both on the blade and in 
the wake. 


That is, by substituting Equation (34-1) in B 
n 
{3 dp = P(nt1) 


where 
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P(3) = a log oa dt By {2b sac) prbe (G5) 
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n-1 
P(nt1) = ae +B, P(n) + E P(n-1) 


where 


iA al 
ie (n-2) a 
B Seen 3) ae 
fo) (n-2) a 
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The integration with respect to ? is performed using the trapezoidal rule. Then the 
integral equation can be represented as a system of linear, simultaneous equations 


for the unknown coefficients a,., b.., and a, 
ij ij it 


' 
1 di K on the blade plane 


PQ 
oy dy itm 5) + D> D>, *vatintoa) Nee 7 (36) 
p=l q=l 


i=1 j=1 k=1 


on the cavity plane 


where the An are the integrations of the double integrals of source distribution in 
Equations (18), (20), (27)-(29) both on the blade and the cavity planes at the &th 
collocation point, and m(ij) is the index corresponding to the coefficient 445° 

By the least-squares method of solving simultaneous equations, a square matrix 


ora is made out of the generally nonsquare matrix Corry where 


B= A. A 37 
eA, 2 


and the right side of the matrix is 


where Co represents the right side of simultaneous Equations (36). The solution of 


the simultaneous equation 


Pe re ety: + S B = B 
D3 pq(ij) 2, &m "pq (km) *« “pq(k) P Oe 
Sd Do LAre 


Pas 
is the desired least squares solution for the cavity source. 


IMAGES FOR THE HUB 


To satisfy the hub boundary condition 


the two image systems for vortex and source distributions were considered separately. 
The hub images for the vortex distribution are the same as those Kerwin’ used. 

That is, the same vortex strength as that distributed on the blade (x,r,,0) is used 

125229) where r, * r, = A Then the trailing vortices tend to 


it 2 H 
cancel the radial velocity on the hub. This is an approximation to a two-dimensional 


at the image point (x 


vortex outside a circular cylinder. 
The hub images for the source distribution need special consideration because 
the cavity sources for a supercavitating propeller are much stronger than those of 
subcavitating propellers. The first approximation of the image source may be taken 
from an approximation of the hub by a sphere. The image system for a point source 
of a unit strength at r = r, of a sphere is a point source at r 


i 
strength ry/ty and a line sink stretching from r = r 


2 % 
ha n/t with the 


9 to the center of the sphere 


with the strength l/r). In this way the total absolute strength of the line sink is 


equal to the point source at r = Lo: Thus, on the sphere, r = ly the normal 
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velocity induced by the source outside of the sphere is canceled by the image source 
system; however, outside the sphere, the velocity field due to the image system 
decays as flow caused by a doublet inside the sphere. 

The image source distribution is multiplied by the same double polynomial as 
Equations (31) and (32), and the hub boundary conditions are considered together with 
the cavity boundary conditions, Equations (28) and (29), to solve the simultaneous 
Equations (33), (34), and (36) by the least-squares method. The solution is checked 
to see if it actually satisfies the hub boundary condition in addition to the cavity 


boundary conditions. 


BLADE SECTION SHAPE 

When the unknown coefficients are obtained from Equation (38), the source 
distribution is computed from Equations (31) and (33). Thus, all the induced veloc- 
ities can be calculated. It is possible to obtain the section shape of the propeller 
blade by integrating the velocity field along the blade reference surface. However, 
this may require a large amount of velocity information. Thus, the method used by 
Renin’ is followed to correct the foil shape derived from supercavitating cascade 
theory in the preliminary design process. At a field point of each section the 


normal velocity component V is obtained on the pressure surface of the blade from 


3n 


the lifting-surface theory, and V is matched with the corresponding normal velocity 


3n 
obtained from supercavitating cascade theory 


m 
ele 
alg eh GP erin GAY (39) 
i oOo 2n 3n 
i=o 
where a, and Cans unknown coefficients 
vee normal velocity component on a foil of a 
supercavitating cascade 
x = distance from the leading edge 


The unknown coefficients are determined by the least-squares method. That is, the 
al : : 
coefficient c_ and the terms a, x are the correction to the normal velocity compo- 
fe) 
nent of the two-dimensional supercavitating cascade caused by the effects of the 


lifting surface, three-dimensional cavity, flow retardation, etc. For a 


Ly 


subcavitating propeller, an airfoil with the same pressure distribution and the same 
F . 3 Gaye ‘ z 

source distribution has been used instead of a supercavitating cascade. When NO 

is represented as the angle of the two-dimensional velocity with respect to the nose- 


tail line and v is represented as the angle of the three-dimensional velocity a 


3n 
with respect to the geometric advance angle 8, then by taking a single term in the 

summation of Equation (39), ay will be interpreted as the corrected angle of attack 
of the nose-tail line with respect to 8, and co will indicate a camber correction of 


the two-dimensional Sodio In the general case the angle of attack will be 


a +—ct a af BOe (40) 


é > : al F F Seah 
and the camber correction due to a; x has to be considered in addition to the 
correction due to co The optimum number of terms in Equation (39) depends upon the 


number of collocation terms. 
FORCES ON THE BLADE 
As in the two-dimensional supercavitating flow, lift and drag should be evalu- 


ated by integrating pressure on the blade surface. The section drag in the direction 


of the nose-tail line is 
cus (pto) diet ioe =o ee (41) 
Dey P V3y Oe 2> SDV DO.) sh DY 
(e) 


since 


Dn dx = 0 (42) 


where Coy is the friction drag, which is the same as Cy for the subcavitating case. 


From the Kutta-Joukowsky theorem, Equation (22), 
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where Vo can now be represented by induced velocities. 


By writing 


where 


Likewise the sectional 
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Thus 
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(43) 


(44) 


Since v, and y are available from the preliminary design calculations 


2n 
iL 
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PL -| x y dx 
oO 
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oe 
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can also be precalculated in the preliminary design. Then the nondimensional thrust 


and power coefficients are obtained in a conventional way: 


Al 
a T Z 
Co = Pas bare i (C, cos -C, sin >) dr 
0) va R ks 
H 
(45) 
i 
_ WQ eZ ; 
Ch a= i ak I r(C, sin +c), cos $) dr 
7 oR We Z 
H 
where T = propeller thrust 
Q = propeller torque 
@ = pitch angle, i.e., 
o6=BR +a 


The efficiency of the propeller is given by 
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NOTE ON PRELIMINARY DESIGN PROGRAM 

As for subcavitating propellers, the preliminary design calculations form a 
separate propeantus however, it is necessary to briefly explain this subject to help 
understand the final design procedures. 

The present preliminary design of a supercavitating propeller is more or less a 
combination of subcavitating lifting-line theory and two-dimensional supercavitating 
cascade enceryso: The cavity drag-lift ratio will be iteratively fed from the latter 
into the former so as to compute the circulation and pitch of the blades needed to 
produce a given thrust. In addition, the foil shape is determined so that the drag- 
lift ratio is not too large for a reasonably thick cavity. Both the cavity thick- 
ness near the leading edge, and the minimum cavity length must be prescribed because 
supercavitating propellers are supposed to have a clean cavity covering the suction 
side of the blade. However, this cavity length is a function of not only foil shape 
but also the cavitation number, which is determined by the design conditions. To 
design such a foil, it is convenient to consider three elementary foils: a basic, 
cambered low-drag foil, such as the two-term camber foil; a flat plate to supply 
angle of attack; and the leading-edge singularity to supply the leading-edge cavity 
thickness. Called a point eae because it produces no lift, the leading-edge 
singularity is especially useful for creating a long cavity. This is so because the 
infinite cavity cavitation number of a supercavitating cascade is linearly propor- 
tional to leading-edge tideinees: see Figure 2. Figure 2 shows the relationship 
between cavitation number and leading-edge thickness for a supercavitating propel- 
le Model 3770 in Figure 2 will be used for the numerical test and is discussed 
in more detail later. 

In general, the specification for the leading-edge cavity thickness is based 
upon the strength requirement. When the cavitation number is relatively large, 
however, the leading-edge cavity thickness may have to be greater than the strength 
condition requires. This greater thickness is needed to allow cavity lengths that 
are more than 50 percent longer than the chord. Each section has different cascade 
parameters, i.e., the solidity and the stagger cece see Figure 3. The present 
preliminary design method allows the option of using either a specified basic load 
distribution or a specified basic camber shape for the blade section. In an actual 


design, the given basic camber shape of the blade section seems to be more conve- 


nient. Given the basic camber shape (the two-term camber in an infinite medium 
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normalized by the lift coefficient, for example) one finds the shock-free angle at 
each station of a supercavitating propeller blade by considering a supercavitating 
cascade of infinite cavity length at each blade section. The angle of attack and 
the point drag are combined with the basic shock-free camber to supply the given 
leading-edge thickness, the minimum cavity length, and the required load distribu- 
tion. In this process there are two options: either the camber or the angle of 
attack per unit lift coefficient is Eee eee a where the former is Case 1 and the 
latter is Case 2. 

The hydrodynamic advance angle B. can be determined by either of two ways: any 
form of tan Bs may be preset, such as tan B. = c/r, or the optimum pitch condition 


for the supercavitating propellers may be preset 


Gewe! lt+e(1-t)- (c+ aA 


dG 


se 
| (46) 
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where € is the drag-lift ratio, and t is the thrust deduction. 

Since the local cavitation number is fixed according to the design conditions, 
and the lift coefficient varies for each iteration, the cavity length may also vary. 
Therefore, the sectional supercavitating cascade problem must be solved for each 
section and for each iteration. This requires a great deal of computer time. We 
resolved this difficulty by conveniently treating the cavity problem as a foil 
having infinite cavity length with a correction for the finite cavity ofrecer The 
former needs to be calculated only once at each blade section, regardless of the 
number of iterations. Only the finite cavity effect is computed for each iteration, 
assuming the same load distribution as for infinite cavity length. Thus, the finite 
cavity effect appears as a reduction of angle of attack and camber. 

The output of the preliminary design program consists of the thrust and power 
coefficients, circulation and lift distribution, tan Bas pitch distribution, 
efficiency, etc. Normal velocity distribution on the cavity and foil, foil cavity 
shape, load distribution, and other necessary data are stored on a tape to be fed 


into the lifting-surface design. 
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NUMERICAL SCHEMES AND COMPUTER PROGRAM 

There are several computer programs for lifting-surface design of subcavitating 
propellers. Supercavitating propellers are similar to wide- and thick-bladed sub- 
cavitating propellers when the cavity thickness is known; therefore, it would seem 
reasonable to use an existing computer program for subcavitating propellers rather 
than to start the entire complicated program from scratch. Among the available 
programs, the recent program by Kerwin was chosen for two reasons: (1) it included 
the effect of sources and vortices, and (2) it could include the effects of rake and 
skew with variable i. 

The main differences between the programs for supercavitating and subcavitating 
propellers are as follows. The strengths of source distributions are not known in 
the supercavitating case because the cavity shape is not known while the blade thick- 
ness is assumed to be known for the subcavitating case. Therefore, the program for 
solving the cavity-source distribution is written according to the present theory 
explained in the previous section and becomes the main frame of the present program. 
Routines in the Kerwin program are used as much as possible. The source is distrib- 
uted not only on the blade but also in the cavity wake according to Equations (31) 
and (34). Therefore, the number of meshes for supercavitating propellers is almost 
as much as two and a half times the number for the subcavitating case. The chord- 
wise load distribution is taken from the supercavitating cascade theory of pre- 
liminary design rather than from airfoil theory. The forces on the blade are ob- 
tained from Equations (41) and (45). 

Figure 4 gives the outline of the flow charts of the lifting-surface program. 

Attention should also be drawn to the following points, which differ from the 
existing program for subcavitating propellers. 

Because of the singular behavior of the kernels of integrals appearing in the 
induced velocity expressions in Equations (17) and (18), all the integrations related 
to bound and trailing vorticity, and the source distribution are integrated analyt- 
ically in a small interval of the radial direction on the blade and cavity where the 
collocation points are located; using Equation (35). 

Because of the singularity along the leading edge, both for source and load 
distribution, the interval that includes the leading edge is treated separately. 
That is, the leading-edge load is the integrated load obtained when the remaining 


load distribution is subtracted from the total load. The leading-edge source 


23 


strength is represented by the normal velocities of two-dimensional cascade results 
solved in the preliminary design program; see Appendix B. 

For convenience, short cavity effects of the supercavitating cascade are com- 
puted in the lifting-surface theory, when the two-dimensional cascade is used as an 
approximate solution to be corrected. 

After computing all the ae coefficients of Equation (37), the simultaneous 
Equations (38) are solved by a subroutine which uses a Gaussian elimination method. 
The accuracy of the solution obtained by the least-squares method is approximately 
checked by a comparison of the left-hand and right-hand sides of Equation (28). Then 
the total source strength on the blade is obtained from Equation (31). The velocity 
component due to sources is obtained from all the velocity component functions 
created for ead in Equation (36). 

Because many existing supercavitating-propeller models have chord lengths 
smaller near the blade tip than near the hub, angular intervals near the tip are too 
small. Thus, the present program allows intervals of all the even degrees, such as 
2, 4, etc. In addition, a l-deg interval is tested for a six-bladed propeller. This 
feature is supplied by using the addition rule of trigonometry, making use of data 
stored for cosine and sine functions. 

When the rake and skew are present in the reference surface, the correct area 
element Hdpdd in all the area integrals has to be taken into account instead of 
H,dpdé. 

The collocation points can be taken to be 10 points on the blade and 5 points 
on the cavity wake. However, if four points are selected on the blade this will be 
exactly the same as the collocation points for the subcavitating case, except for 
the extra points on the cavity wake. 

In addition to other preliminary data, the ae and Qa in Equation (44) are pre- 
pared in the preliminary design and are conveyed through an input tape. As in the 
subcavitating case, friction is taken into account, but only on the pressure side. 

The output routine is created in a format similar to that of Kerwin by con- 
sidering data taken from the supercavitating cascade of the preliminary design for 
the blade-section and the cavity shapes. 

The input to the computer program includes the same information as is used for 
the subcavitating propeller design, such as rpm; ship speed; propeller diameter; hub 


diameter; helical distance from an arbitrarily fixed reference plane, distance from 
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the x-y plane to the leading and trailing edges of the blade; inflow velocity ratio 
V,/Vo3 number of collocation points on the blade, cavity wake, etc. Also, as 
mentioned previously, the data tape created by the preliminary design program for 
all the necessary sectional data from the supercavitating cascade theory is fed into 
the lifting-surface design. 

The output consists of the components of induced velocities resulting from the 
vortex and source distributions, correction factors for source strength and camber, 
angle of attack relative to 8, pitch distribution, thrust and torque coefficients, 


and efficiency. 


NUMERICAL EXPERIMENTS FOR SOLUTION 

The cavity model, the solution for cavity source strength, the number of inter- 
vals for vortex and source distributions, the number of collocation points, and the 
like cannot be determined theoretically; instead, they must be determined according 
to the behavior of numerical output. The output, especially, should show conver- 
gence, which could be built into the program, if the program is simple and does not 
require too much time. However, a large program, such as the present one, which 
requires considerable computer time, cannot be run for all values of parameters so 
as to check convergence of the solution for each design. Instead, it may be enough 
to check several aspects of parametric changes for a typical case and to assume that 
the other cases will reasonably follow the typical case. 

In the present work, designing the Model 4717 supercavitating propeller, various 
convergence checks have been performed, and the results will be shown in the follow- 
up. The design conditions for Model 4717C are shown in Table l. 

The cavity truncation locations were varied from 1.5 to 2.8 chord lengths to 
check the cavity model. Figures 5 and 6 show that the truncation at 2.2 chord 
lengths and at 2.5 chord lengths produces almost the same pitch and camber distri- 
butions (differing less than 1%); the pitch distribution is about 3% larger than at 
1.8 chord lengths. 

The influence of the choice of angular interval on pitch and camber distri- 
bution is also shown in Figures 5 and 6. By changing 2-deg intervals to 1-deg 
intervals, the camber correction factor se increased about 6% and the pitch diameter 
ratio P/D increased about 1%. Since the magnitude of camber is so small, a change 


in the correction factor of 10% is within the manufacturing error. If the blade tip 


25 


has more than five intervals, 2-deg intervals are sufficient; finer intervals, which 
increase costs significantly, appear to be unnecessary. 

To check whether there were enough collocation points, we chose 10 points and 
created the foil shape by plotting streamlines. The results have been compared with 
the case of four collocation points given in Figures 7 and 8. The simple approxi- 
mation obtained using Equation (39) is shown to express amazingly close agreement 
with the plotted streamlines. 

In Equation (28) Vo is the local mean speed, which is not known without exami- 
nation, so the value from the lifting line theory has been substituted for it. How- 
ever, a more accurate formulation would be Equation (53) (Appendix A) instead of 


Equation (28). 


u u V 
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The computer results for these two cases were almost the same. The present program 
used Equation (47), although slightly more computer time was needed. When G/V. is 
quite large, it may improve the solution. 

To check whether the solution satisfies the boundary conditions well, the left- 
hand and right-hand sides of Equation (47) were plotted in Figures 9-12; and the 
radial components of velocities on the blades were plotted, Figures 13 and 14. These 
calculations were made with and without the hub boundary conditions being satisfied 
when the degrees of p and x in the double polynomials in Equations (31) through (33) 
were taken as 3 in one case and as 5 in the other case. During this process, we 
noticed an interesting phenomenon: an instability occurred in the numerical value 
of radial velocity for the solution which did not satisfy the hub boundary condition. 
That is, if the hub boundary condition was not specified, a slight change of parame- 
ters, such as cavity length, number of intervals, or the degree of polynomials, 
produced large changes in radial velocities. Yet, the numerical values of the thrust 
and torque coefficients or the pitch distribution did not change too much. This may 
be because the linear boundary conditions on a cavity or foil do not include any 


constraint on the radial velocity. In the present problem, the only constraint on 
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the radial velocity is on the hub condition. Thus, the hub boundary condition is 
needed not only to find the hub effect on the pitch distribution but also to make 
the solution stable. 

The specified boundary conditions are well satisfied in general, although when 
the boundary condition on the hub is included, the cavity conditions are slightly 
less accurate, as shown in Figures 11 and 12. The differences in the radial veloc- 
ities occurring for the case with and the case without the hub boundary condition 
(see Figures 13 and 14) indicate an instability. The radial velocity satisfying the 
hub boundary condition in Figure 13 is the stable solution. In Figure 15 the pitch- 
diameter ratio is shown for Model 4717C with hub images. In Figure 16 the camber 
correction factors are shown for Model 4717C with and without hub images. When the 
computed results obtained without consideration of the hub boundary condition 
happened to have radial velocities with small values, the results were very close to 
the solutions obtained when the hub boundary condition was considered. The numerical 
results reported in the following discussion were obtained without satisfying the 
hub boundary condition. In the cases given, the instability phenomenon was not 


noticed. 


NUMERICAL EXAMPLES FOR PROPELLER DESIGN AND DISCUSSIONS 

Many supercavitating propellers have been designed and tested in the past; of 
that number, two propellers, DINSRDC Models 3770 and 3870 were hosen- to determine 
if this design program is reasonable. The former propeller has three blades and a 
low advance coefficient, and the latter has four blades and a high advance coeffi- 
cient. Experiments showed that both propellers had smooth cavities. The experi- 
mental results and the previous design calculations are available. 

The design and performance characteristics of the two propellers are shown in 
Table 2. 

It is extremely difficult to compare the present numerical results to the 
experimental results for propellers that were designed using an entirely different 
method. The present program is intended for design, not prediction. The present 
program does not produce data on leading-edge cavity thicknesses, input that is 
essential to the design of propellers similar to Models 3770 and 3870. To check 
the reasonableness of the present program we guessed at the leading-edge cavity 


thicknesses for Models 3770 and 3870; this is presented in Figure 17. Because the 
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actual cavity thickness was never measured, this is not a scientific estimate. The 
leading-edge cavity thickness selected for Model 3770 results in an almost infinite 
cavity length at every section of the blade, except near the tip and the hub. How- 
ever, because the design cavitation number of Model 38/70 is not very small, the 
selected leading-edge cavity thickness is not large enough to induce a smooth sheet 
cavity all over the blade. Therefore, the leading-edge cavity thickness of Model 
3870 was corrected to give a cavity length at least 50 percent longer than the chord. 
To do this, an extra leading-edge point drag was added to the cascade theory used in 
the present design method, as explained previously. The design thrusts used in the 
present lifting-line computations are the experimental values listed in Table 2, 
where Case 1 is for cases without angle of attack and Case 2, with preset angle of 
attack. It is seen that the efficiencies of the propellers, as predicted by the 
present method, are very close to the measured efficiencies, even in the preliminary 
design stage of calculations. 

Pitch distributions obtained from the preliminary design and lifting-surface 
design computations, according to the two design approaches, Case 1 (without angle 
of attack) and Case 2 (with angle of attack) are shown in Figures 18 through 21 
together with the pitches of the two propeller models. The pitch values obtained 
from the preliminary design calculations are higher than those obtained from lifting- 
surface calculations because, in preliminary design, the effect of flow retardation 
is not considered. The pitch distribution is also related to the leading-edge 
cavity thickness. In general, when the leading-edge thickness increases, the pitch 
also increases; however, the efficiency decreases slightly. The pitch distributions 
for the predictions and for the models are noticeably different. This is because 
the optimum lift distribution and the pitch angle, which are influenced by the blade- 
cavity interference in the present method are quite different from those of the 
models. If these factors are taken into account, all results appear reasonable 
compared with those of the models. 

The lifting-surface corrections to the source distribution for the two pro- 
pellers are shown in Figure 22. Although the correction for Model 3770 is close to 
1, it is about 20 percent greater near the trailing edge than at the leading edge for 
Model 3870. If the correction factor is unity, this means that the cascade source 


strength is the same as the blade cavity source strength. This source strength has 
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the main influence on pitch, thrust, and efficiency. The calculated efficiency is 
close to that obtained in the model experiments, as shown in Table 2. 

The lifting-surface camber correction factors for the two propellers. are shown 
in Figures 23 and 24. The camber correction factor for Model 3870 is much larger 
than for Model 3770, as in the calculation of Venning and Bab ermeanis The pitch 
distribution and camber correction curves for Model 3870 are quite different from 
those for Model 3770. The former has shorter cavities with four blades and large 
expanded area ratio (EAR), and the latter has longer cavities with three blades and 


smaller EAR. 


EXPERIMENTAL EVALUATION OF THEORY 

Finally, the Center conducted an experimental program to evaluate the present 
method for designing supercavitating propellers. Two supercavitating propellers 
were designed, using the present me chad re for a 200-ton (181 metric ton) hydrofoil 
craft. The propeller design criteria are given in Table 3. The propeller design 
characteristics are given in Table 2. 

Two model propellers were manufactured from these designs. The geometry of 
these propellers is given in Table 4 and drawings of the propellers are shown in 
Figures 25 and 26. 

The experimental program was divided into two phases. The results of the first 
phase, the measurements of blade-cavity shapes, are reported in the following 
section, while the results of model propeller performance are described in the last 


section. 


EXPERIMENTAL INVESTIGATION OF BLADE-CAVITY 
THICKNESS DISTRIBUTION 


A series of experiments was performed to determine how well linear theory 
predicted the upper cavity surface location for Propellers 4717C and 4738A. For 
these experiments, brass pins of varying lengths were attached to the backs of the 
propeller blades. During propeller operation in the 36-in. variable-pressure water 
tunnel, one could see when the pins came into contact with the upper cavity surface. 

This experimental procedure has already been used to verify the upper cavity 


19520 


surface location for Propeller 4699, and the parent design of Propellers 4738A 
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and 4717C. In addition, cavity heights have been measured in similar ways for a 


supercavitating flat plate and a Tulin-two-term section by Christopher and Jonmeonke 


DESCRIPTION AND LOCATION OF PINS 

Number-four brass machine screws with heads cut off were used as pins. These 
pins were attached to the blade backs by drilling and tapping holes perpendicular to 
the surface at 12 locations. The upper cavity surface location was defined as a 
point on a line, perpendicular to the nose tail line, that runs through the center 
of the tapped hole at the blade surface; see Figure 27. The holes were drilled and 
tapped perpendicular to the back of the blade to cant the pins slightly away from the 
reference line (the line perpendicular to the nose tail line). A slight error was 
introduced, but the machining process was greatly simplified. The locations of these 
pins were as follows: 10, 30, 60, and 90 percent of chord at nondimensional radii 
(r/R) of 0.361, 0.544, and 0.726; see Figure 28. Since the blade at 10 percent of 
chord was too thin to tap, the pins at these locations were soldered in place. The 
pins at all other locations were screwed in and secured by tiny electrical lock nuts 
(see photos in Reference 19). 

Three sets of pins were used for testing Propeller 4717C. When installed, the 
first set of pins protruded above the back of the blades to a height that corre- 
sponded to three times the distance from the back of the blades to the theoretically 
predicted upper cavity surface. The second set of pins protruded by a factor of 1.67 
and the third set by a factor of 1.0, the latter being the theoretically predicted 
cavity height. However, the pins at 10 percent of chord varied from this order; 
their heights corresponded to cavity-height factors of 3, 2.33, and 1.0. The pins at 
the 10 percent of chord locations had to be filed by hand to the correct height, and 
the factor of 2.33 rather than 1.67 was used to ensure that the pin-height would 
exceed the experimental cavity thickness. The experiments have shown, however, that 
the theory overpredicted the cavity height near the leading edge, and the pins could 
have been filed to a height corresponding to a factor of 1.6/7. 

Four sets of pins were used for Propeller 4738A. The first set of pins pro- 
truded above the back of the blade to a height corresponding to a factor of 1.8 times 
the distance from the back of the blade to the theoretically predicted cavity sur- 


face. The other three sets corresponded to factor of 1.4, 1.0 (theoretical) and 
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0.6. The pins at the 10 percent of chord locations were filed to a height corre- 
sponding to the multiplication factor used for the other pins in the set being 
tested. 

A thin coat of international yellow paint was applied to the tips of the pins 
prior to testing to aid in visual observation. The experimenters could thus locate 


the pins easily when the propeller was revolving. 


EXPERIMENTAL PROCEDURE 

The propeller rpm was increased until the cavity enclosed all the pins. This 
procedure was begun with the set of brass pins that protruded highest. While the 
pressure and velocity in the 36-in. variable pressure water tunnel was held constant, 
corresponding to the design cavitation number, o = 0.34, the propeller rpm was 
gradually reduced in decrements of 10. Each time the propeller rpm was reduced, a 
hand-held strobe unit was used to observe visually whether the pins were in or out 
of the cavity. These observations could be made rapidly, although the pins that 
barely touched the upper cavity surface required more attention than did the others. 
When all pins protruded through the cavity surface, the procedure was repeated with 
another set of pins. A depth micrometer was used before and after each test to 
measure the height of the pins above the back of each blade. This was done to ensure 
that the pins were at the correct height and had not moved during testing. 

As the pins began to break through the cavity surface, a furrow or small 
groove formed in the surface, accompanied by some spray or cavitation behind the 
pin. A directional strobe unit with variable light intensity made these furrows 
much more visible. This phenomenon has been recorded in several color photographs. 

Large, international-yellow numbers painted on the backs of the blades proved 
invaluable during testing. Also, each propeller hub was coded with a series of dots, 


the number of which corresponded to the number on each blade. 


DISCUSSION OF EXPERIMENTAL RESULTS 

Figures 29-31 compare linear theory predictions of cavity height with experi- 
mental results for Propeller 4717C. These figures show three experimental upper 
cavity surfaces corresponding to three values of J, one of which is the design J 


(1.037). Note that the following relationship gives the advance angle, 8, at each 


radial blade section, 


3: 


B= Cane (V,/21rn) 


cares Guo 


From this relationship, one can determine, approximately, the corresponding shifts 
in upper cavity surface that result from small changes in angle of attack. For 
example, in Figure 29, at the design J (1.037), 8 is determined to be 42.44 deg. 
For a J value of 1.0, B = 41.40 deg. Therefore, a change of J corresponding to 
0.037 has caused, approximately, a one-deg change in angle of attack, which shifts 
the cavity surface upwards as shown. 

Theoretically, the section lift and cavity thickness for Propeller 4717C are 
generated entirely by camber and point drag (note the blunt nose in Figures 29 
through 31). That is, no incidence was used in the design to generate lift or cavity 
thickness. In Figures 29 though 31, the theoretical prediction of cavity height 
agrees fairly well with the experimental data. Near the leading edge, however, the 
theory appears to overpredict cavity thickness. Also, visual observations indicated 
that the backs of the blades at all radial sections on Propeller 4717C were wetted 
to about 2- or 3-percent of chord from the leading edge. At this point, separation 
was caused by a locally flat area that was inadvertently machined onto the back of 
the blade. Although this local flat was almost microscopic, it effectively caused 
separation. Apparently, very near the leading edge some portion of the blade metal 
was interfering with the upper cavity streamline. 

Figures 32 to 34 compare linear theory predictions of cavity height with experi- 
mental results for Propeller 4738A. Note in these figures the large amount of point 
drag or blunt nose indicated by the theory. This results because both Models 4738A 
and 4717C were designed to have approximately the same full-scale stress levels. 
This dictated that the maximum, theoretical, cavity thicknesses for Models 4738A and 
4717C would be almost the eamencs To obtain the same maximum thickness in a shorter 
distance, we used a large amount of point drag together with incidence and camber to 
generate the theoretical cavity. 

Figures 32 to 34 show three experimental upper cavity surfaces corresponding to 


three values of J. Note that at r/R = 0.361, Figure 32, the blade was fully wetted 


at the design value of J (1.037); therefore, cavity heights for three other values 
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of J have been shown. At a J value of 0.98, the experimental cavity surface coin- 
cides with the predicted cavity surface height at aft locations on the blade. This 

J value represents an approximate increase in angle of attack of 1.6 deg over the 
design value of J. For the cavity surface, corresponding to a J value of 0.96, the 
increased incidence is equal to about 2.2 deg. However, as mentioned before, the 
cavity did not spring from the leading edge. It moved down the span as rpm was in- 
creased. Since the water-tunnel velocity was held close to 35 fps (10.688 m/sec) 

and since the model propeller diameter was 16 in. (40.64 cm), the difference in model 
rpm corresponding to the J values of 0.98 and 1.037 was 88. This corresponds to an 
increase of about 58 in full-scale rpm. It is also interesting to note that, accord- 
ing to performance evaluation experiments, an increase of 29 rpm, over the 1000 rpm 
of full-scale design would give the design thrust. 

At the two outer radial positions, r/R = 0.544 and 0.726, full cavitation did 
occur at the design J (1.037), but the theory overpredicted the cavity surface 
height. As with Propeller 4717C, the back of the blade near the leading edge of 
Propeller 4738A was wetted to about 2- or 3-percent of chord. 

To understand more fully the discrepancy between theory and experiment, the 
reader should recall that a point drag is a linear theoretical model of the leading 
edge cavity thickness represented by a point singularity. Experimental results indi- 
cate that the actual separation point at the leading edge must be carefully chosen, 
for example, as a slope discontinuity of the blade surface, to achieve the designed 
leading edge cavity thickness; if the predicted leading edge cavity, not just 70 
percent of it, had been filled by a material up to 2 percent of chord from the 
leading edge, the experimental results would have almost coincided with the theory, 


except very near the hub, where the hub effect is important. 


CAVITATION PERFORMANCE CHARACTERISTICS OF SUPERCAVITATING 
PROPELLERS 4717B, 4717C AND 4738A 


BACKGROUND 

Propeller 4717C was originally manufactured as Propeller 4717B. Propeller 4717B 
was identical to Propeller 4717C except for the backs of the blades, which had a 
shape to conform to the predicted cavity shape at design operating conditions. The 


primary purpose of Propeller 4717B was to determine, by observation, how well the 
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blade section shape represented the blade cavity shape. Following characterization — 
and observation, Propeller 4717B was finish cut to the final design version, Pro- 
peller 4717C. 

Propeller 4738A is a six bladed propeller with the same expanded area as the 
previous four bladed propeller, and is designed for the same conditions as the other 


propellers. 


EXPERIMENTAL PROCEDURE 

Cavitation performance characteristics and cavitation observations were ob- 
tained in the 36-in. variable pressure water tunnel. Tunnel water velocities were 
measured by the tunnel venturi system. The scope of the experiments is iver in 
Table 5. 

Tunnel pressure and water velocity were set to establish each cavitation number 
and then propeller revolution rate was varied to cover a range of advance coeffi- 
cients. Propeller thrust and torque were measured at each condition and sketches 
were made of the cavitation present. The Reynolds number, R_, during the experiments 


5 6 iy 
ranged from’'7.5 X10} ‘tol 5.6): *' 10%. 


PRESENTATION OF DATA AND DISCUSSION 
The thrust and torque data were reduced to nondimensional coefficients of Ki and 
K,. Propeller efficiencies were calculated from faired values of K,, and K.. The 


lene tee performance characteristics of the three propellers are ce in 
Tables 6 through 8. 

Curves representing the faired data, from Tables 6 through 8, are shown as an 
example in Figure 35. Curves of maximum-speed thrust loading (K,/I°) have been 
added to the performance curves for Propeller 4738A. The intersection of the ng 
curve and the Ko curves at the design sigma (0) determines the predicted operational 
point for each propeller. A comparison between the design operational points and 
the points predicted by the experimental data is given in Table 9. 

Sketches of the back cavitation present on the propellers at two cavitation 
numbers are given in Figures 36 through 38. These sketches cover a range of advance 
coefficients from partially cavitating to fully cavitating conditions. If advance 
coefficients had lower values than those shown, the propellers, at the same cavi- 


tation number, would also be fully cavitating. With only one exception, propellers 
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contained no face cavitation over the range of cavitation numbers and advance co- 
efficients covered. The sole exception was at an advance coefficient of 1.2, where 
some leading edge face cavitation was observed at cavitation numbers of 0.75 and 


lower. 


CONCLUSIONS 

At design speed coefficient and design 0, Propeller 4717B contains practically 
no back cavitation. If advance coefficient is reduced slightly, at design o, the 
backs of the blades are covered by sheet cavitation from the blade tip to 50 percent 
radius. This indicates that the predicted cavity shape over this part of the pro- 
peller blades is quite accurate. 

Neither of the designed propellers, 4717C and 4738A, had face cavitation at the 
design operational points. Propeller 4717C essentially had full back cavitation and 
Propeller 4738A had back cavitation from about 35 percent radius to the tip of the 
blades at the design operational point. 

The propeller theory slightly overpredicted the available thrust for both pro- 
pellers. Propeller 4717C would require 6 percent more rpm and 8 percent more power 
than predicted to reach design speed. Propeller 4738A would require 5 percent more 
rpm but 4 percent less power than predicted to reach design speed. If the operating 
point is defined as the speed and rpm where the propellers absorb the available 
maximum power, Propeller 4717C would operate at Vi = 58.7 knots and rpm = 1016, and 
Propeller 4738A would operate at ae = 61 knots and rpm = 1054. The propeller 
efficiencies at these conditions are 66 percent and 6/7 percent, respectively. 

It has been neeemaie dae that the nonlinear effects on lift and drag of cavi- 
tating foils are approximately equal to -0.5 c?/ (140) and -0.5 Cyc, / (+0). There- 
fore, propellers designed according to the linear theory would produce less thrust, 
as indicated in the experimental results. However, the leading edge cavity thickness 
was slightly smaller than the design thickness due to the unmatched separation point 
and 70 percent filling of cavity thickness. Thus the drag should have been a little 
less than predicted by the linear theory, as was found experimentally for Propeller 
4738A. Because of the decreased efficiency of Propeller 4717C, however, more care 
may be required in calculating the cavity drag due to the blunt leading edge. 


All in all, the design theory predicted the cavity thickness and the propeller 


performance quite well at the design point, within the bounds of error to be 
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expected of linear theory. It should be remembered that the present design theory 
does not include semiempirical formulae or factors. In addition, this verifies that 
the point drag of linear theory is very useful in solving one of the important 
problems relating to supercavitating propellers--how to make propeller blades having 


sufficiently thick leading edges without paying too much of an efficiency penalty. 
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Figure 1 - Coordinate Systems for Blade Surface Reference 
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Figure 2 - Example of Relation Between Choked Flow 
Cavitation Number and Leading-Edge Thickness of 
Supercavitating Cascade for Blade Sections of 
Supercavitating Propeller Model 3770 
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Figure 3 - Example of Cascade Parameters for 
Supercavitating Propeller Model 3770 
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READ INPUT: PROPELLER GEOMETRY, V,/V, 
RPM, ETC., FROM CARDS 

READ SUPERCAVITATING CASCADE DATA OF 
PRELIMINARY DESIGN FROM A TAPE 
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Figure 4 - Flow Chart for Lifting-Surface Design of 
Supercavitating Propellers 


39 


6 = ANGULAR INTERVAL 
2= CAVITY LENGTH 


(5,2) = (2, 1.8) 


(5 = 1 deg, & = 1.8) 


“0 0.2 0.4 0.6 0.8 1.0 
r/R 


Figure 5 - Computed Camper Correction Factors co of Supercavitating 


Propeller Model 4717 for Different Angular Intervals 6 and 
Cavity Truncation Points 2% 
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Figure 6 - Computed Pitch-Diameter Ratio of Supercavitating Propeller 
Model 4717 for Different Angular Intervals 6 and Cavity Truncation 
Points & 
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Figure 7 - Computed Shape of Blade Face of Model 4717 at r/R = O. 
for Two Numbers of Collocation Points 
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Figure 8 - Computed Shape of Blade Face of Model 4717 at r/R = 0.3978 


for Two Numbers of Collocation Point 
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Figure 9 - Cavity Boundary Conditions without Hub Image and with Cubic 
Polynomial Source of Model 4717 
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Figure 10 - Cavity Boundary Condition with 5-Degree Polynomial and 
without Hub Image of Model 4717 
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Figure 11 - Cavity Boundary Condition with Hub Images and Cubic 
Polynomial Source of Model 4717 
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Figure 12 - Boundary Condition on Cavity with 5-Degree Polynomial and 
Hub Images of Model 4717 
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Figure 13 - Radial Component of Velocity with Cubic Polynomial 
Source of Model 4717 
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Figure 14 - Radial Component of Velocity with 5-Degree Polynomial 
Source without Hub Image of Model 4717 
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PRELIMINARY DESIGN 


FINAL DESIGN 
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Figure 15 - Computed Pitch-Diameter Ratio of Model 4717 
with Hub Image 
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WITHOUT HUB-IMAGE 
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Figure 16 - Camber Correction Factor cy of Model 4717 with 


Hub Image, 6 = 2 Degrees and & = 2.5 
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MODEL 3870 


MODEL 3770 
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Figure 17 - Design Cavity Thickness t/(cC, ) 
at Leading Edge x/c = 0.1 
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LIFTING-LINE DESIGN 


MODEL 3770 
LIFTING-SURFACE DESIGN 
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Figure 18 - Pitch-Diameter Ratio of Model 3770 with 
Leading-Edge Point Drag (Case 1) 
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LIFTING-SURFACE DESIGN 
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Figure 19 - Pitch-Diameter Ratio of Model 3770 without 
Leading-Edge Point Drag (Case 2) 
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Figure 20 - Pitch-Diameter Ratio of Model 3870 (Case 1) 
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Figure 21 - Pitch-Diameter Ratio of Model 3870 (Case 2) 
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Figure 22 - Lifting-Surface Source Correction Factors (Case 2) 
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Figure 23 - Camber Correction Factor co of Model 3770 
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Figure 24 - Camber Correction Factor co of Model 3870 
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Figure 27 - Placement of Brass Pin Perpendicular to Back of Blade 


58 


LEADING EDGE 
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LEADING EDGE 


Figure 28 - Pin Locations at 10, 30, 60 and 90 Percent of Chord at Non- 
Dimensional Radii (r/R) Values of 0.361, 0.544, and 0.726 for 
Propellers 4717C and 4738A 


59 


OL 


OLTLy APeTTedorg A0; 
T9€°O = U/A We sadeys ARTAeD TeoTAtdwy pue [TeoTZeACeY, Jo uostzedwoy - 67 eAn3Ty 


(3/X) GHOHD 4O NOILIVYS 
60 80 £0 9°0 g0 v0 £0 c0 LO 


(AHOSHL) ALIAWS YaddN 
30V18 4O NOV 


GLL=*°V7 


SLNIOd VLVG TwiNawivadxa { 4€0t =F C) 
OL=rO 


0:0 


0° 


OL 


00° 


oc 


so'0- 


0 


0 


0 


0 


(3/A) LHDISH ALIAWS 


60 


OL 


6 


0 


80 £0 


Javi1g 40 nova isco 


OLTLy AeTTedorg 103 
'4G¢°O = U/a 2e sedeys AqTAeD TeoTAtTdwy pue [eoTJeAOeY], Jo uosTieduoD - QE sANsTYy 


(9/X) GHOHD 4O NOILOVYS 


9°0 


G0 


v0 £0 20 LO 00 
g0'0- 


1 00'0 


S0'0 


OL'O 


(AYOSHL) ALIAWS YAddN 


SLO 


zoL=r O 
SLNIOd VLVG IWLNAWIHadxa { ZeoL =r 
soL=r V7 02'0 


(9/A) LHDISH ALIAWS 


61 


OL 


OLTLy AeTTedorg 103 
97L°O = U/A 3e sedeys AjTAeD TeoTAtTdwY pue TeOTIeA0s8y], FO uostieduoj -— [E eAnsTy 


(9/X) GHOHD JO NOILOVYS 
60 80 £0 9°0 s'0 v0 £0 20 LO 


30Vv1¢8 40 NOV 


oLl=r 7 


SLNIOd VLVG 1WLNAWIY3dxS zoL=f O 
zeot=r 


00 


(9/A) LHDISH ALIAWS 


62 


OL 


V8ELy AeTTedoig 105 
T9E°O = U/a 2e Sodeys AQTAeD) Teotatdwuy pue TeoTJeA0sYI 


(3/X) GHOHD 40 NOILOVUYS 
60 80 £0 90 ‘0 v0 


LE0'L =f :‘NOILIGNOD NDISAG 
Lv LAM 4GV14 JO WOVE 


jo uostieduoy - 7E sanstTy 


£0 c'0 LO 


(AYOSHL) ALIAWS Yaddn 


oeo=" 7 
SLNIOd VLVG TVLNAWIHAdXS { 96:0=F O 
s6o=r TC) 


0o'0 


(9/A) LHDISH ALIAWS 


63 


OL 


V8ELY teTTedoig 103 
7G°O = U/at 7e sodeys Aqtae) Teotatdwy pue TeorjJez109y] Jo uostiedwoy) - ¢¢ eANnsTy 


(9/X) GHOHD 4O NOILOVYS 
60 80 £0 9°0 G0 v0 £0 c0 LO 0'0 


S0°0 


Javid 40 WOVE 


(AYOSHL) ALIAWS Y3ddN 


OL="° V7 
SLNIOd VLVG TWLNaWIH3dxd { ZOL=FO 
feo ler) 0z'0 


(9/A) LHDIAH ALIAWS 


64 


OL 


V8Ely AeTTedorg 105 
97Z°O = U/A Je Sadeysg ARTAeD TeoTATdwy pue TeoTJerOey], Jo uostTieduog —- HE VsANnSTy 


(3/X) GHOHS 40 NOILOVYSA 
60 30 £0 90 G0 v0 £0 c0 LO 


Aqvia 4O 1ova 


(AYOSHL) ALIAWS YaddN 


OL="° 7 
SLNIOd VLVG IVLNAWIYadXa } ZOL=F O 
zeor=e O 


00 
g0'0- 


00°0 


s0'0 


OL'0 


SLO 


0c'0 


(9/A) LHDISH ALIAWO 


65 


cL 


io 


V8ELy AeTTedoig Jo ooueWIOFI0g UOTIeATAeD - GE eANsTy 


f LN3I91I44309 GAONVAGV 
LL OL 60 80 £0 9°0 G0 


£80°0 = 2f/+» re 


‘0 


ee 


0 
GL 
= SS = = 


SL0 


= 
oO 


N 
ro) 


OY 
So 


Ss 
fo) 


oe 
r—) 


oS 
3 


(O> OL) LN319144309 ANDYOL GNv (+) 1N319144909 LSNYHL 


66 


o = 0.34 


Figure 36 - Sketches of Cavitation Present on the Back of Propeller 4717B 
at Two Cavitation Numbers 
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J=0.9 1.0 1.04 1.1 


Figure 37 - Sketches of Cavitation Present on the Back of Propeller 4717C 
at Two Cavitation Numbers 
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J=0.9 1.0 1.04 1.1 


0 = 0.34 


Figure 38 - Sketches of Cavitation Present on the Back of Propeller 4738A 
at Two Cavitation Numbers 
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TABLE 1 - PROPELLER DESIGN PREDICTIONS 


Propeller Number 4717C | 4738A ; 


58.84 knots (30.35 m/sec) 
1000 
0.092 


IG (0)S}7/ 


6 
1.410 
0.34 
0.67 
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TABLE 2 - DESIGN AND PERFORMANCE CHARACTERISTICS OF 
SUPERCAVITATING PROPELLERS 


Propeller 


Z 
P/D (0.7) 
EAR 
oi (Os 7) 
J 


Experiment 
K 


ae 
n 
by Venning and Habermas 
Preliminary Deaton 


(Lifting-Line Theory) 
Lifting-Surface Design 


Case 1** 
Case 2*% 
by Venning and Hebeemena 
Preliminary Deeden 
Case 1** 
Case 2** 
Lifting-Surface Design 
Case 1 


Case 2 


“Team and Haberman. 


*A corrected value of Ke as given in 
Reference 2. 


*kCase 1 - without angle of attack; Case 
preset angle of attack. 
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TABLE 3 —- PROPELLER DESIGN CRITERIA 


Power, maximum continuous 


1-t 
1l-w 
t 
Shaft Centerline at 
Design Speed 


Maximum D 

No. of Blades 

Design Objective 
Minimum Rotative Speed 


Hump Thrust Margin 
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16,000 hp (11,931 kW) 
for two propellers 


0.925 
0.870 


6.82 ft (2.079 m) 


5 ft (1.524 m) 
4, 6 
maximum speed 


750 rpm 


20 percent 


TABLE 4 — MODEL PROPELLER GEOMETRY* 


Propeller Number 


4717¢ 4738A ) 


Design Parameters 


Z 4 6 
Diameter, in. (cm) 16.000 (40.64) 16.000 (40.64) 
Pitch,** in. (cm) 22.656 (57.546) 22.560 (57.302) 


P/D** 1.416 1.410 
EAR 0.495 0.492 
Chord length,** in. (cm) 4,560 (11.582) 3.017 (7.663) 


*Interpreted from drafting room offsets. 


**At O.7R. 


TABLE 5 - SCOPE OF EXPERIMENTS 


*In feet per second. 
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TABLE 6 — CAVITATION PERFORMANCE OF PROPELLER 4717B 


SIGMA (o)* = 3.000 SIGMA* = 1.500 


*Cavitation number. 
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TABLE 7 - CAVITATION PERFORMANCE OF PROPELLER 4717C 


SIGMA* = 3.000 SIGMA* = 1.500 


“t 


*Cavitation number. 
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TABLE 8 - CAVITATION PERFORMANCE OF PROPELLER 4738A 


SIGMA* = 6.000 SIGMA* = 3.000 
1OKQ mK 


SIGMA* = 1.500 SIGMA* = 0.750 


- 5000 
-5500 
-6000 
-6500 
- 7000 
- 7500 
-8000 
-8500 
-9000 
-9500 
1.0000 
1.0500 
1.1000 
1.1500 
1.2000 


SIGMA* = 0.500 SIGMA* = 0.340 


-6000 
-6500 
- 7000 
7500 
-8000 
- 8500 
-9000 
-9500 
1.0000 
1.0500 
1.1000 
1.1500 
1.2000 


*Cavitation number. 
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TABLE 9 — PROPELLER OPERATING POINTS 


Propeller Number 
4717B ATT ei Waren 


Performance 
Parameters 


*Rounded from 58.84. 
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APPENDIX A 
CAVITY BOUNDARY CONDITION 


From the Bernoulli equation with respect to coordinates fixed on the blade 


ae] 


2 


ao v 
+ 5) = + Om (48) 


ol 


ae 


where P and q are the pressure and the speed of a fluid particle relative to the 
blade. Thus 


P P ae ac 
(ga gece 
fe) fe) 2 2 + F= pe 
(49) 
=G: Vo 
Therefore 
V 
QL 
mew ve oy 
7 P Ve s s 
Again from the Bernoulli equation 
RP 2 
Gy OW iL Zine, 
a oh s 5 ((V+u,) tu} 
or 
u= ae 
cf 
fey 2VG ela) NG (IT at “7 
poet tee Aa oD 
spas s s Ss 


This is Equation (25) or (26) 
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Writing 


Wij, = AW se Cea ep oy. (52) 
Equation (28) can be written 
ue be V 
BSE ACA) ( Pepe Rema OMe iat 
V ZEN, Vv Vv 2V Vv 
Ss s s 
or 
m 
u V u 2 u Vv 
Techs hy eC ase ae fea gis 
(+e 72\--§ (1+ sf) - Vv TOy (53) 
s s s 
m 
G aes m 
if v is as small as y ory terms including G ur and G u, can be neglected. 


s 
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APPENDIX B 
LEADING-EDGE SOURCE 


We assume the normal velocity on the foil cavity near the leading edge has the 


form 


We solve for ays ay» and a, with three points of x 


3) 
xuVOe O25), O05) and) Oe 
a. —— + a. + a. 0.025 = b 
it V0.025 - 3 z 
a Su 8 a, + a, 0.05 = by 
V¥0.05 
ili 
‘a ap aia. ola (gdb = Io) 
Then with the determinant 
ONS 
V0.025 
me = A O40s 1 O,0s0R706 
V0.05 
Me ne ee 
V0.1 


we obtain 
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ie) 
WW 


(0.05b -0.075b,+0.025b,)/D 


at 1 
ay = (-0. 2891b) +0. 5533985b,-0.2044243b,)/D 
eh (1. 3098582b)-3.1622776b,+1.8524193b,)/D 


The total flux near the cavity is 


x 


Wey) dx = mx 


oO 


mx is the flux of the corresponding wedge (or uniformly distributed Ve up to x). 


That is 
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